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ABSTRACT 

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function 
established? Such intermediate and high-mass stars may be born from relatively massive pre-stellar gas 
cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes 
such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their 
surrounding clump medium. Their internal pressure is provided by a combination of turbulence and 
magnetic fields. On the other hand, the Competitive Accretion model requires strongly sub- virial initial 
tIh I conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate 

Qh, and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we 

have identified four prime examples of massive (~ 100 Mq) clumps from mid-infrared (MIR) extinction 
mapping of Infrared Dark Clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H + in 
these objects, which are consistent with them being starless. Here we present ALMA Cycle observations 
of these 4 clumps that probe the N2D + (3-2) line at 2.3" resolution. We find 6 N2D + cores and determine 
their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the 
predictions of the Turbulent Core model of virialized, magnetized (with Alfven Mach number rriA ~ 1) 
and self-gravitating cores that are bounded by the high pressures of their surrounding clumps. However, 
pT) ■ in the most massive cores, with masses up to ~ 60 Mq, our results suggest that moderately enhanced 

t^J- ' magnetic fields (so that vtla — 0.3) may be needed for the structures to be in virial and pressure 

equilibrium. Magnetically regulated core formation may thus be important in controlling the formation 
of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar initial mass 
function. 

■ Subject headings: ISM: clouds, dust, extinction — stars: formation 
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1. INTRODUCTION 

The two main theories for massive (and intermediate-mass) star formation, Core Accretion (e.g. McLaughlin & Pudritz 
^ ; 1996; McKee & Tan 2003, hereafter MT03) and Competitive Accretion (e.g. Bonnell et al. 2001; Wang et al. 2010), invoke 
Uh . very different initial conditions for the gas about to collapse to form a high-mass star. The Turbulent Core Accretion 
_ model of MT03 starts with massive near-virial-equilibrium starless cores that will collapse to form individual stars or close 
binaries. These cores can be considered to be scaled-up versions of the much more common low-mass cores known to form 
low- mass stars (Shu et al. 1987). The main differences from low- mass cores are that the pressure support must be provided 
by nonthermal forms (turbulence and/or magnetic fields) and the typical environments where massive cores and stars 
form, self-gravitating massive clumps (i.e., proto-star-clusters) , are at much higher pressure. Competitive Accretion also 
involves fragmentation of massive gas clumps, but now into protostellar seeds with initial masses only of order the thermal 
Jeans mass — typically much less than a solar mass under these high pressure conditions. Those protostellar seeds that 
happen to be in high density regions then later accrete previously unbound gas to eventually become intermediate-mass 
and high-mass stars. The numerical simulations of Bonnell et al. (2001), from which the competitive accretion theory was 
developed, involved global, near free-fall collapse of a gas clump, which either started at or evolved to a very sub-virial 
state. Krumholz et al. (2005) showed such a sub-virial state was needed for the competitive accretion rate to be large 
enough to be relevant for massive star formation, i.e. allowing formation within timescales < 1 Myr. Thus to distinguish 
between these theoretical models we need to find and measure the dynamical state of intermediate-mass and high-mass 
starless cores: how close are they to virial equilibrium? 

Answering this question is challenging because massive starless cores are rare, typically far away, small in angular size, 
usually surrounded by large quantities of other molecular gas, and expected to suffer depletion of many molecular species, 
such as CO and CS, that are often used to measure the mass and kinematics of molecular clouds. Furthermore, it is 
difficult to infer the degree of support a core receives from large scale magnetic fields. 
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1.1. Target selection of intermediate-mass and high-mass starless cores 

To find massive starless cores we started by studying infrared dark clouds (IRDCs) (e.g. Perault et al. 1996; Egan 
et al. 1998) using Spitzer IRAC 8/zm GLIMPSE survey (Benjamin et al. 2003) images. We developed a mid-infrared 
extinction (MIREX) mapping technique to derive mass surface densities, E, of clouds, probing up to E ~ 0.5 g cm -2 (i.e., 
iV H = 2.1 x 10 23 cm" 2 ) and on scales down to 2" (Butler & Tan 2009, hereafter BT09; Butler & Tan 2012, hereafter 
BT12). The method depends on the dust opacity but not the dust temperature. Figure[TJi shows E maps of 4 core/clumps 
(CI, Fl, F2, G2) selected from the larger sample of 42 studied by BT12. The properties of these objects are also listed 
in Table [TJ Selection of these particular sources was guided by them having large values of E, still being dark at 24 and 
70 /im (Spitzer MIPSGAL (Carey et al. 2009) images were analyzed), and being in relatively quiescent environments with 
respect to other star formation activity (MIR sources are absent in a ~ 20" diameter aperture around the E peak). 

1.2. CO Freeze-out and Deuterium Fractionation 

In these cold (T ~ 15 K; Pillai et al. 2006), high density (jih > 10 5 cm" 3 ; BT12) environments, freeze-out of molecules 
such as CO on to dust grains is expected, as has been observed in low-mass starless cores (e.g. Caselli et al. 1999). 
Hernandez et al. (2011) found widespread CO depletion in the filamentary IRDC H of the BT12 sample. It has also been 
reported in other similar clouds by, e.g. Zhang et al. (2009) and Fontani et al. (2012). Molecules suffering freeze-out are 
poor tracers of starless core kinematics. However, as CO is depleted from the gas phase, the level of deuteration of certain 
species is expected to increase (Dalgarno & Lepp 1984). This is because CO reacts with and thus destroys H2D + (that 
has formed by reaction of Hj with HD), so when CO is depleted the concentration of H2D + builds up and so deuteration 
of other species still present in the gas phase, such as N2H + , becomes enhanced. 

In fact the 4 IRDC clumps we selected were found to have the highest levels of deuteration fraction, Z?f rac = N NoD + /N-n 2 h+ 
among the 10 massive starless clouds studied by Fontani et al. (2011), with D frac = 0.38, 0.43, 0.40, 0.70 for CI, Fl, F2, 
G2, respectively. This compares with a mean Df rac = 0.12 for the other 6 sources in the sample, which were not selected 
from MIR extinction. This indicates that the MIREX-based selection method is efficient at finding regions that are cold 
and dense, and not yet forming stars. Note that the protostellar sources in the sample of Fontani et al. tend to have 
significantly lower values of Df rac — 0.04. 

Deuterated molecules have been shown to be the best kinematic probes of the coldest and densest conditions associated 
with low-mass star formation (e.g. Caselli et al. 2002; Crapsi et al. 2007). We will apply the same methods in the more 
extreme environments of IRDCs, where there tends to be a larger and more pressurized mass of cold, dense gas, compared 
to low-mass star-forming regions such as Taurus. We will thus use N2D + to both identify pre-stellar cores — efficiently 
distinguishing them from their surrounding clump material — and then as their kinematic tracer. 

1.3. Predictions of the Turbulent Core Accretion Model of Virialized Cores 

MT03 modeled clumps and massive pre-stellar cores as singular polytropic virialized spheres in hydrostatic equilibrium, 
with power law density, p oc r -fep , and pressure, P oc r~ fep , distributions with fiducial value k p — 1.5 (for both clumps 
and cores), implying kp = 2(k p — 1) = 1.0. Such values are consistent with those derived from MIREX mapping by 
BT12, who found fc PiC i ~ 1.1 (for clumps) and k PtC ~ 1.6 (for cores). MT03 assumed cores have a boundary with surface 
pressure, P s , c , that is in approximate equilibrium with the pressure of the immediately surrounding environment of the 
star-forming clump. This pressure was assumed to be related to the mean pressure in the clump, P c \, via P s c — </>_p iC P c i, 
where <f>p. c is a constant of order unity. The value of (j)p,c depends on the location of the core within the clump: MT03 
assumed cores were typically at radial locations of 0.3i? c i so that <pp. c ~ 2 (for A; PiC i = 1.5). In this paper we will define 
the clump properties from the material that is more localized around each core: we set R c \ = 2R C and evaluate its mass 



Table 1 

IRDC Core/Clumps Observing Targets 



Core name a 


RA(J2000) 


Dec(J2000) 




d c 


Dt b 

J - / irac 




(h m s) 


(° ' ") 


(km s - 1 ) 


(kpc) 




G028-C1 


18:42:46.9 


-04:04:08 


+78.3 


5.0 


0.38 


G034-F1 


18:53:19.1 


+01:26:53 


+57.7 


3.7 d 


0.43 


G034-F2 


18:53:16.5 


+01:26:10 


+57.7 


3.7 d 


0.40 


G034-G2 


18:56:50.0 


+01:23:08 


+43.6 


2.9 


0.70 



a From BT12 

b From Fontani et al. (2011) 

c Kinematic distance from Simon et al. (2006) 

d Kurayama et al. (2011) report a parallax distance of 1.56 kpc for 
IRDC F, although the reliability of this has been questioned by Foster 
et al. (2012). 
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surface density, £ c i, as the average value (observed via the MIREX maps) in the annulus from R c to 2R C . In this case, 
4>p,c = (1 — kp /3)(i? c / 'R c i)~ kp —> 4/3 in the fiducial case of k p = 1.5, kp = 1 (and 1.07 in the case of k p = 1.1, kp = 0.2). 

MT03 assumed the mean pressure in the clump was set by its self-gravitating weight, i.e. P c \ = (frpGYf^, where 4>p — 
(37r/2O)/g0geom<^,Bavir,ci, where f g is the gas mass fraction of the clump, gC om is a numerical factor of order unity that 
accounts for the effect of nonspherical geometry, 4>p accounts for the effect of magnetic fields and a V i r ,ci = 5(cr^ 1 )i? c i/(GM c i) 
is the virial parameter of the clump. MT03 adopted a fiducial value of f g = 2 /3 (with the remainder composed of a nascent 
stellar cluster), but for the IRDC clumps we consider here, we will assume f g ~ 1. Aspect ratios of 2:1 lead to </> gC0 m — 1-1, 
so such elongations are a relatively minor effect and, like MT03, we assume ^ ge0 m = 1- For magnetic field strengths such 
that the Alfven Mach number via = 1, MT03 showed that (f)p = 1-3 + 1.5m^ 2 — > 2.8. Below, for our fiducial estimates, 
we will consider the possibility of a range 0.5 < vtia < 2, corresponding to 7.3 > <fip > 1.7. Finally, we assume that the 
clump is virialized with a V i r ,ci = 1- With these values cf>p — 1.32 (c.f. 0.88 in MT03). 

The radius of a core in virial equilibrium, including pressure equilibrium with its surroundings, is (MT03) 

A V"/ M, ' ^- ->- ,/2 



R "" = 0071 Kw*^;) Iti^J » c w 

where A = (3 — k p )(k p — l)f g — > 3/4. The mass average velocity dispersion of a virialized core is related to that at the 

surface via t7 c ,vir = [2(3 — k p )/(8 — 3fc p )]<7 CiS — > (6/7)er CiS . Note that a c — c c /0^ 2 , where c c = {P c / p c ) 1 / 2 is the effective 
sound speed in the core, so that 

2(3 -k p ) ( M£ \ 1/8 / M c \ 1/4 / S cl \ 1/4 



Mr. 



1/4 / v \ 1/4 



1.09 — ^ J kms" 1 . (4) 

V 60A/ m y V i g cm- 2 y V ' 

o 

2.355a v i r ) versus size relation: 



gem 

The above expressions can also be combined to give a velocity dispersion (or FWHM line width, AiVir = (81n2) 1 / 2 cr 



V 1 1 



,2(3 - k p ) f <t> P ,c<l>p\ 1,A ( E d \ 1/2 fRo ■ xV2 



2.27 v ■ L1 9 ™i k -i (5) 

1/2 X 1/2 



^-►1.44 — ^ kms' 1 . (6) 

\1 g cm 1 ) \0.1 pc / 

We will now compare these predictions, in particular for <7 c , v ir and i? c ,vir, with our sample of massive starless cores. 

2. OBSERVATIONS 

We used ALMA in Cycle compact configuration to observe N2D + (3-2) (231 GHz, Band 6, 2.3" angular resolution, 
~ 9" maximum detectable scale, 27" field of view, 0.08 kms -1 velocity resolution) with a single pointing for each core, 
centered at the peak E of the BT12 MIREX map, and sharing one track amongst the 4 sources with a total integration time 
of about 1.0 hour per source. Several other species, DCO + (3-2), DCN(3-2), 13 CS(5-4) and SiO(5-4) were also searched 
for in the same spectral setup. From the line-free regions of the four spectral windows we measured the continuum flux 
with an average wavelength of 1.338 mm. The values of the RMS noise per beam in CI, Fl, F2, G2 were 0.272, 0.213, 
0.336, 0.174 mJy/beam, respectively. 



3. RESULTS 

3.1. Identification of starless cores via N2D + emission 

The integrated intensity maps of N 2 D+(3-2) and DCO+(3-2) of the clumps CI, Fl, F2, G2 are shown in Fig. [Tb & 
c. Strong detections of both lines are seen in all cases. DCN(3-2), 13 CS(5-4) and SiO(5-4) were not detected, even from 
stacked spectra towards the N2D + cores. Continuum emission from the line-free regions of the four spectral windows, 
with mean wavelength of 1.338 mm, was detected and is shown in Fig. Id. 

The critical density of N 2 D+(3-2) is 2.9 x 10 6 cm" 3 , while that of DCO+(3-2) is 3.5 x 10 6 cm- 3 . However, the emission 
from DCO + tends to be more widespread. We expect this is because DCO + (formed via CO + H 2 D+) is following the 
more extended distribution of CO, which is likely to suffer significant depletion in the cold, high-density cores, where 
N2D + is enhanced. There is generally good correspondence between the morphology of the line emission and structures 
seen in the MIREX maps. For example, the "V" -shaped structure of DCO + emission in Fl is also seen in the MIREX 
map. Several, even quite weak DCO + features in G2, also show themselves via MIR extinction. 

Some, but not all, localized MIREX column density structures reveal themselves via N2D + (3-2). Following the results 
of studies of low-mass starless cores described in JTJ we define these to be starless core candidates, finding 6 objects 
(Cl-N, Cl-S, Fl, F2, G2-N, G2-S) that have emission exceeding the 3a noise level within I — b — u-space. Positional 
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CI (N&S) Fl F2 G2(N&S) 




2fl".329 2B°.327 2B°.32S 28°.323 28°,321 34*.424 34°.422 34*42 34°.41B 34°.4I6 M°-44 34*.43° 34°.43B 34°.434 34°.432 34°.785 34°.783 34°.7fl1 34°.77e 34°.777 



Fig. 1. — The four massive starless core/clumps observed by ALMA: columns from left to right show CI, Fl, F2, G2. (a) First row: MIREX 
S maps in gem -2 (BT12). The inner circles show the deconvolved extent of the N2D + cores identified in row (b), while the outer circle has 
a radius twice as large and shows the region used to estimate the surrounding clump envelope, (b ) Second row: ALMA Cycle observations 
of N2D + (3-2) integrated intensity. Six cores, Cl-N, Cl-S, Fl, F2, G2-N & G2-S, are defined by their 3<r contours in position- velocity space 
Note that not all high S regions show strong N2D+emission, but N2D + cores do correlate well with structures seen in the MIREX maps, (c) 
Third row: ALMA observations of DCO + (3-2) integrated intensity, which is generally more widespread than N2D + . (d) Fourth row: ALMA 
observations of the 1.34 mm dust continuum emission. 

core boundaries are defined using the projection of this 3cr contour, after deconvolving with a Gaussian equivalent to the 
ALMA synthesized beam. The positions and average radial sizes, i? c , are listed in Table [2] The diameters of the cores 
are all larger than the angular resolution of the observations. Core Cl-N appears to exhibit some substructure in its 
N2D + (3-2) emission, whereas the other sources appear to be relatively monolithic. 

The morphology of the continuum emission also generally matches that seen in the molecular line emission, especially 
that traced by N 2 D+(3-2). 

3.2. Core kinematics and velocity dispersion 

The integrated N2D + (3-2) spectra of the cores are shown in Fig. [2] These spectra are a blend of 47 hyperfme components. 
We fit model spectra to these data to derive the centroid velocity, Vlsrn 2 d+; anc ^ the observed ID velocity dispersion, 
°n 2 d+ obs; a l so listed in Table HJ This modeling allows for the possibility of optically thick parts of the hyperfme line 
complex, but we find all spectra can be well-modeled assuming optically thin line emission. We assume a gas temperature 
of T = 10 ± 3 K (see § 13.3. 2\i to remove the thermal component of the line to thus assess the nonthermal component of 
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Fig. 2.— N 2 D+ (3-2) spectra (y LSR -frame) of the six identified cores: Cl-N, Cl-S, Fl, F2, G2-N & G2-S, shown at the maximum resolution 
of 0.08 kms -1 . The red continuous line shows the best fit model profiles. 



the velocity dispersion. We then sum this in quadrature with the sound speed of the molecular cloud, assuming the same 
temperature value of 10 ± 3 K and a mean particle mass of /x = 2.33to p , to derive the total core velocity dispersion as 
derived from N 2 D + (3-2) emission, on 2 d+- 

We do not use DCO + to measure core kinematics since it is expected to be somewhat under abundant in CO-depleted 
regions and thus likely to preferentially trace the extended core envelopes. Nevertheless we present in Table[5]the values of 
Vlsr,dco+ an d CTE>co+,obsj m ^ e g ra t e d over the same region of the N2D + -defined core. The centroid velocities of N2D + (3-2) 
and DCO + (3-2) are very similar (always within 0.2 km s" 1 ), while the observed velocity dispersion of DCO + tends to be 
larger by up to ~ 35%, which is probably related to this species tracing a larger region around the core. 

In Fig. [3]we show the first moment maps of the 6 cores as traced by N2D + (3-2) and DCO + (3-2). We integrate over a 
velocity range of ±2.5 kms -1 centered on the mean velocity of the N 2 D + cores Cl-S, Fl, F2, and G2-S. In the CI region, 
cores Cl-N and S show a velocity offset of 1.8 kms -1 , but with some N 2 D + emission at intermediate velocities. In G2, 
the cores are within 0.4 kms -1 of each other. On the larger scales of the DCO + (3-2) emission the first moment maps 
generally reveal quite disordered structure, with multiple discrete features spanning a velocity range of a few km s _1 . The 
exception is Fl, which shows quite coherent velocity structure across its "V" -shaped morphology. 

3.3. Core and clump mass surface density, mass and density estimates 

3.3.1. Estimates from MIREX maps 

We use the BT12 MIREX maps to estimate the mass surface density of the clump material, S c i, immediately surrounding 
the identified cores, from annuli extending from R c to 2R C . This choice of size of the annulus follows the method of BT12 
(see also ijl.31) . We find values of £ c i ~ 0.2 — 0.3 g cm -2 (Table [2]). We assume these estimates have a 30% uncertainty 
due to the choice of opacity per unit gas mass in the IRDC material (BT12 assume K$nm — 7.5 cm 2 g -1 based on the 
moderately coagulated thin ice mantle model of Ossenkopf & Henning 1994). However, in very high column density 
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CI (N&S) Fl F2 G2(N&5) 




Fig. 3. — The four massive starless core/clumps observed by ALMA: columns from left to right show CI, Fl, F2, G2. (a) First row: 
ALMA Cycle observations of the first moment map of N2D+(3-2) emission, integrated over a 5 km/s velocity range centered on the mean 
N2D+core velocities of Cl-S, Fl, F2 and G2-S. (b) Second row: ALMA observations of the first moment map of DCO+(3-2) emission in the 
same velocity reference frames as in row (a). 



regions, the MIREX method of BT12 will begin to underestimate the true value of S, when the background intensity 
that makes it through the cloud becomes comparable to the instrumental noise in the image. BT12 refer to these as 
"saturated" regions. Only F2 is formally classified as "saturated" by BT12, but this is based on a global definition within 
the IRDC and does not allow for modest variations in the foreground intensity across the cloud. 

Butler, Tan & Kainulainen (BTK), in prep., have examined archival Spitzer- IRAC data on IRDC C that probe to 
significantly smaller instrumental noise levels than the Spitzer-GLIMPSF, images analyzed by BT12. We have evaluated 
E c i from these BTK maps and find values that are 35% and 25% larger for Cl-S and Cl-N, respectively. We will perform 
our dynamical analysis for the following cases: (1) uniformly using the BT12 maps for all 6 cores; (2) replacing the BT12 
maps with the BTK maps for cores Cl-S and Cl-N (these entries are shown in square brackets in Table [5]). 

We then use the MIREX maps to evaluate the maximum mass surface density of the core, £ c ,max> averaged inside 
R c . In one sense this is an upper limit for the core properties because it assumes all the line of sight material in the 
IRDC is associated with the core. However, because of potential saturation in the MIREX maps, our estimates of the 
total mass surface density through the IRDCs in these regions could be underestimates. From the BT12 maps, we find 
S c ,max — 0.2 — 0.4 g cm" 2 , only slightly higher than the values of S c i. This suggests that either only a small fraction of 
the total column of material along the line of sight is associated with the core (with the rest being part of the surrounding 
clump) or that the mass surface densities of the cores are underestimated because of saturation in the MIREX maps. To 
assess this latter possibility, for the CI cores we also evaluate £ c ,max from the BTK map, finding values that are 56% and 
47% larger for Cl-S and Cl-N, respectively, than those derived from the BT12 maps. In the BTK map, these cores have 
values of £ c ,max that are 21% and 23% larger than E c i (compared to 8.5% and 13% in the BT12 maps). As a result of 
the uncertainty in assessing what fraction of £ c ,max is associated with the cores, after first deriving core properties based 
simply on £ c ,max, we will also consider two further methods to estimate E c . 

Utilizing S c max , the "maximum" core mass is derived as M C:max = E Cimax 7ri? 2 . We assume 20% distance uncertainties 
(we have adopted the kinematic distance estimates of Simon et al. 2006), which leads to 50% total uncertainties in M Cimax 
(summing errors in quadrature). However, IRDC F has a reported astrometric distance of 1.56 kpc (Kurayama et al. 
2011), only 42% of the kinematic distance. We thus perform a separate analysis for the F cores using this closer distance 
(these entries are shown in rounded brackets in Tabled]). Note however that Foster et al. (2012) have called into question 
the validity of the astrometric distance, since the data quality on which it is based are relatively poor and the implied 
streaming motions if it is correct are very large (~ 30 km s _1 ). 

The derived maximum core masses, M Cimax , are in the range from ~ 1 to ~ 60 Mq. Cl-N & S have similar masses at 
the upper end of this range (using the BTK maps) and have the potential to form massive (i.e. > 8Mq) stars. Next, G2-N 
has about 9 Mq, followed by the F cores (assuming the kinematic distance) and then G2-S. These may be the progenitors 
of intermediate-mass or even low-mass stars. Note, however, that these star-forming environments are at significantly 
higher mass surface densities, and thus pressures, than most previously studied regions of low-mass star formation, such 
as in the Taurus molecular cloud. 

With £ c ,max and R c we can also estimate the maximum volume density assuming spherical geometry. Table[2]lists values 
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of the maximum total H nuclei number density as riH,c,max = 3M Ciinax /(47ri? 3 /XH), which given the above assumptions 
have 36% uncertainties. These values range from a few xlO 5 cm -3 to few xlO 6 cm~ 3 . 

We now attempt to account for the overlying clump material that is not associated with the cores. First, we assume this 
material has the same column density as the surrounding clump evaluated from R c to 2R C . We expect that subtracting 
off this clump envelope leads to a mimimum estimate of the core mass, since if the local IRDC geometry is sheet-like 
or filamentary this would lead to an overestimate of the clump material that is along the line of sight to the core. We 
evaluate £ c ,min = £ c ,tot — £ c l and list its values in Table [2j We also then derive M Cim ; n and nH, c ,min- This method of 
envelope subtraction substantially reduces the estimated core mass by factors of about five to ten. The most massive 
cores, Cl-S & N now have masses of about 10 Mq (using the BTK maps). These estimated densities are < 10 5 cm~ 3 . 

3.3.2. Estimates from 1.34 mm dust continuum emission 

Additionally, we utilize the observed 1.338 mm continuum emission to estimate core properties. Our ALMA Cycle 
compact configuration observations are sensitive to angular scales from ~ 2" to 9", so will tend to pick out core rather 
than clump material. The total mass surface density corresponding to a given specific intensity of mm continuum emission 
is 



E mm = 6.03 x 10- 3 [j£L_ j { o.OlcrnH- 1 ) 1338 ^ i L075T ^^ 33 «J " X J gCm < ^ 

where A1.338 = A/1. 338mm and T^io = T^/IOK. At this frequency, our preferred choice of k„ = 5.95 x 10 _3 cm 2 g _1 , based 
on the moderately coagulated thin ice mantle dust model of OH94 and assuming a total (gas plus dust)-to-refractory- 
componcnt-dust-mass ratio of 141 (Draine 2011). OH94 dust models that are more coagulated have k„ that is about 23% 
larger. Overall we adopt a 30% uncertainty in k„. 

An estimate of the dust temperature is also needed. From observations of ammonia inversion transitions, Pillai et al. 
(2006) measured IRDC temperatures ~ 15 K. However, we expect such observations tend to probe the lower density 
envelopes around the IRDC cores, as the inversion transitions of NH3 have critical densities of only ~ 10 4 cm~ 3 . The dust 
temperature can also be constrained from the FIR spectral energy distribution. For example, fitting Herschel PACS and 
SPIRE data, Peretto et al. (2010) found temperatures of ~ 10 — 12 K in the central region of an IRDC. 

Here we use the fact that some IRDC cores are seen in absorption at wavelengths as long as ~ 100/im to place constraints 
on the dust temperature. This method has the advantage of allowing us to derive constraints on scales of the Herschel 
angular resolution, i.e. down to 6" for the 70/xm band. Approximating the cloud as a ID slab that has a total optical 
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Fig. 5. — Dependence of i^i/T^o with E in the Herschel-P ACS 70, 100, & 160 (im wavebands for IRDCs near I = 30° and for various dust 
temperatures (the dotted lines show IK increments between the labeled temperatures). OH94 moderately coagulated thin ice mantle dust 
opacities have been adopted. Note the cloud needs to be cold (< 15, 12, 9 K) to appear dark at 70, 100, 160 fim, respectively. 



depth tv at frequency v and integrating the radiative transfer equation yields 

Iu,i = Iu fi e- T » + f " e-^-<\j v /K v )drl, (8) 



n 



where I Vi \ is the transmitted intensity towards the observer from the near side of the cloud, I v> q is the background 
incident intensity on the far side of the cloud, j u is the emissivity of the cloud material (i.e. its dust), and k v is the cloud 
opacity. Approximating the dust in the IRDC as being in thermal equilibrium and at a constant temperature so that 
j v /n v B v {T d ) = (2hp 3 /c?)(e hu / kTd - 1)"\ the Planck function, we have 

^i =e -- + ^) (1 _ e - n (9) 

The background intensity can be estimated empirically from the observed emission in the Galactic plane. Li & Draine 
(2001) considered Galactic infrared emission in the MIRS region (44° < I < 44°40', -0°40' < b < 0°). Integrating 
their model over the Herschel-PACS bands at 70, 100, 160 /im gives intensities I v = 500, 1270, 1690 MJy sr^ 1 at mean 
wavelengths of 74.0, 103.6, 161.6 ^m, respectively. At I — 30°, b — 0°, Bernard et al. (2010) estimated "offset" intensities 
of 688.0 and 1982.3 MJy sr _1 for the PACS 70 and 160 /xm bands, which provides an approximate estimate for the 
intensity of the diffuse emission in this direction. Using an average scaling factor of 1.27, the 100 /im flux at I = 30° 
is estimated to be 1620 MJysr" 1 . Our target IRDCs C, F, G are at I = 28.3°, 34.4°, 34.8°, so we will use the results 
for I ~ 30°. BT09 estimated that for IRDCs C, F, G the fraction of diffuse emission that is in the foreground of the 
cloud is 0.266, 0.193, 0.14, respectively. Adopting a value of 0.2, so that 80% of the observed intensity is the value of the 
background intensity behind the IRDCs, we have I v p = 550, 1290, 1590 MJy sr _1 . 

With the above fiducial values of Figure [S] shows as a function of mass surface density, E = t„/k„, for 

different dust temperatures, T^. We have evaluated I v ,x/Iv,q for the fluxes that would be received in the 70, 100, 160 /im, 
i.e. the Herschel-PACS wavebands, integrating over the filter response function, the Li & Draine (2001) spectrum of the 
diffuse Galactic background emission, and over the OH94 opacity function (we find mean opacities of At7o^ m = 1.14cm 2 g _1 , 
KiooAim = 0.603 cm 2 g _1 and Kieo^m — 0.290 cm 2 g _1 for the OH94 moderately coagulated thin ice mantle dust model). 

We examined images of the IRDCs in the Herschel data archive. IRDC C has been imaged at 70, 100 & 160 /im. At 
70 ijluy, CI appears globally dark (i.e. relative to the low-intensity diffuse emission beyond the BT09 ellipse, not just 
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locally dark with respect to its immediate surroundings). At 100 /jm it has a similar, perhaps slightly lower, intensity 
as the faintest parts of the IRDC surroundings, while at 160 /im it appears moderately brighter (although is still locally 
dark). Comparison with Fig. [5] suggests a (mass-weighted line of sight) temperature of ~ 10 K at CI. IRDC F has images 
at 70 & 160 /im. Fl & F2 are locally dark at 70 fim. While F2 also appears to be globally dark, this does not appear to 
be the case for Fl due to the relative proximity of bright sources. At 160 fim, the intensities towards both Fl & F2 are 
brighter than the diffuse surroundings. IRDC G has images at 70 & 160 /im. G2 appears globally dark at 70 /j,m, but 
brighter than the IRDC surroundings at 160 /tm. 

Note that we expect the temperature in the N2D + cores to be lower than the above constraints from ~ 100 /j,m 
shadowing, which are mass-weighted averages along the line of sight through the IRDC. Studies of lower-mass starless 
cores have seen temperatures down to ^ 6 K in the very central regions (e.g. Crapsi et al. 2007; Pagani et al. 2009). 
Given the above constraints, we adopt a dust temperature of Td — 10 ± 3 K for the N2D+ cores. As we will see, the 
densities are high enough that we expect the gas and dust temperatures to be well-coupled, so we set T = T<j. 

We then evaluate S c mm via eq. [7] Since the uncertainties introduced by temperature are quite large and asymmetric, 
we show upper and lower uncertainty bounds separately (we have combined uncertainties in quadrature for the upper and 
lower sides separately). The values of E c ,mm are always greater than E C)ln j n . They are less than S Cjmax , except in the case 
of Cl-S, where it is about 20% greater than the BTK E c ,max value. These values still agree given the uncertainties. Also 
recall that £ derived from MIREX mapping of saturated cores is a lower limit, and this may well be the case with Cl-S. 

We then proceed to derive masses (M Cjlnm ) and densities (nH, c ,mm) based on E mm . We will tend to regard the core 
masses derived from mm dust emission as the most accurate measure, rather than M c>max or M c , m in- This is because the 
intcrferometric observations of the mm continuum emission filter out contributions from structures > 9", which are the 
scales expected of the surrounding clump. The cores themselves have diameters ranging from 3.7" to 7.3". Note, that 
-Wc,mm still has uncertainties of about a factor of 2, due to assumed temperature (±3 K), opacity (30%) and distance 
(20%) uncertainties. 

3.4. Dynamical state of the cores 

For each estimate of core mass, i.e. M c max , M Cim - ln and M Cjlnm , we calculate a c y n - from equation iQ}. In each case, 
we first consider the properties of the cores and their clump envelopes as derived from the BT12 MIREX maps. We also 
show the results for Cl-N and S when using the BTK extinction map, and the results for Fl and F2 when adopting the 
near distance of 1.56 kpc (see Table [J). In addition to the uncertainties in er c , v j r due to errors in M c and £ c i, we also 
allow for a range of magnetic field strengths in the core, such that 0.5 < itia < 2, i.e. 7.3 > 4>b > 1-7 (see £|1.3p . We then 
evaluate the ratio of the observed velocity dispersion derived from the N2D + (3-2) spectrum to the prediction from virial 
equilibrium, <Tn 2 d+ /cc,vir- We average these values for the sample of 6 cores (assuming the uncertainties are uncorrelated, 
although in reality there are likely to be correlated systematic uncertainties, e.g. affecting mass determinations due to 
choice of k or X^). We present three averages: first with core properties derived from the BT12 maps, second with the 
CI cores evaluated using the BTK maps (and the remaining cores with the BT12 maps), and third being the same as the 
second case, but with the F cores evaluated at the near distance. We regard the second case as being the most accurate. 
We repeat the above analysis for i? c .vir (evaluated from eq. (2) and the ratio R c /R c ,vir (note the distance uncertainty does 
not affect R c / R c , v - n ). 

Quoting results for the second case, for M C)lnax the core sample has 0N 3D +/er C)V i r — 0.670 ± 0.086, i.e. moderately 
sub-virial, while R c /R c y lr ~ 1.03 ±0.17. Cores would appear to be sub-virial if their masses have been overestimated. For 
-Wc,min) the sample has cr N2D + /cr c , vir — 1.21 ±0.16, i.e. moderately super-virial, and i? c /i? CiV ir — 3.28 ±0.54. Cores would 
appear to have sizes larger than their predicted equilibrium sizes if their masses have been underestimated. Finally, for 
M Cimm the sample has <7n 2 d+ /°c,vir — 0.813q;^3 and R c /R c ,vii — l-54};26- Given the estimated uncertainties, we conclude 
that the sample of cores, with masses estimated from mm continuum emission, have properties broadly consistent with 
the virial equilibrium predictions of the Turbulent Core Model. There is a tentative evidence that the cores have velocity 
dispersions that are marginally sub-virial, but as we shall see this is degenerate with the assumed magnetic field strength 
in the cores. The core sizes appear to be slightly larger than the equilibrium size. We note that this result may be 
influenced by our, somewhat arbitrary, choice of defining the core size via the 3a N2D + (3-2) contour. 

Focusing on the most massive core, Cl-S, with M c , mm ~ 62.52§ 9 8 M©, we find that <7N 2 D+/°'c,vir,mm — 0.477o;3|f and 
Re/ Rc. vir. mm ~ 0.897J 

!522- This may indicate that the core is quite strongly sub-virial, and so should be undergoing 
fairly rapid global collapse (its free-fall time is only ~ 50, 000 years - Table [3]). Alternatively, the core could be closer to 
virial equilibrium and supported by stronger large-scale magnetic fields than are assumed in the fiducial Turbulent Core 
Accretion model. 

The assumed fiducial value of the Alfven Mach number rriA = V3a c /vA = V3cr c / {B / ^Ai:p c ) = 1 implies a mean 
background field strength in a core of 

B c = yi^^ml 1 = 297ml 1 (^^) ^ (g^pi) A«G. (10) 

In Table [3] we summarize the dynamical properties of the cores (based on the mm continuum core masses, the BTK 
extinction map for Cl-N & S, and the far distance for IRDC F), including B c for itia — 1, evaluated using eq. [TU] with 
o~ c = <?n 2 d+- These values are ~ 100 — 300 pG. 

To assess what field strength would be required for the core to be in virial equilibrium with <7n 2 d+ /o"c,vir,mm — 1) note 

that 

^c,vir,mm 4* b (this assumes that the same value of cipplics in the clump envelope, which influences <^p ; els in 
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Table 2 
Identified N 2 D + Cores 



Core property (% error) 


Cl-N 


Cl-S 


Fl 


F2 


G2-N 


G2-S 


Average 


I (°) 


28.32503 


28.32190 


34.41923 


34.43521 


34.78097 


34.77838 




b (°) 


0.06724 


0.06745 


0.24598 


0.24149 


-0.56808 


-0.56829 




C (") 


o.oo 


O.Ul 


2 49 


1 87 

-L .O / 


^ (\7 
0.\) i 


2 31 




d (kpc) (20%) 


5.0 


5.0 


3.7 
(1.56) b 


3.7 
(1.56) b 


2.9 


2.9 




R c (0.01 pc) (20%) 


8.18 


8.75 


4.46 
(1.88) 


3.35 
(1.41) 


5.16 


3.25 




Ken N n+ (kms -1 ) 


81.18±0.03 


79.40±0.01 


56.12±0.01 


57.66±0.04 


41.45±0.01 


41.80±0.01 




°"N 2 D+,ob s (kms- 1 ) 


0.367±0.032 


0.365±0.016 


0.172±0.011 


0.376±0.047 


0.288±0.013 


0.240±0.015 




°"N 2 D+,nt (kms -1 ) 


U.oOoztU.Uo.Z 


n q«i -Un ni « 


n 1 fi/i _i_n m o 


n Q7o_un c\a i 
V.o ( Z±U.U4 i 


U.zooztU.Ulo 


n OQ/i -Un ni ^ 




CT N 2 D+ (kms- 1 ) 


0.409±0.031 


0.407±0.019 


0.250±0.021 


0.417±0.044 


0.340±0.018 


0.300±0.020 




V LSR,DCO+ (kms^ 1 ) 


81.18±0.02 


79.59±0.02 


56.26±0.01 


57.53±0.02 


41.52±0.01 


41.91±0.02 




^DCO+.obs (kms- 1 ) 


0.452±0.021 


0.422±0.019 


0.233±0.004 


0.513±0.017 


0.285±0.012 


0.266±0.016 




S rl (fcm^) (30%) 


0.340 


0.289 


0.217 


0.324 


0.214 


0.194 




[0.459] a 


[0.360] a 












Sc.max (gem- 2 ) (30%) 


0.369 
[0.554] 


0.326 
[0.442] 


0.248 


0.355 


0.225 


0.212 




M c , max (M ) (50%) 


37.1 
[55.8] 


37.5 
[50.9] 


7.41 

(1.32) 


5.97 
(1.06) 


9.01 


3.36 




"H.c.max (10 b c m - d ) (36%) 


4.68 
[7.03] 


3.87 
[5.24] 


5.78 
(13.7) 


11.0 
(26.1) 


4.53 


6.78 




°"c,vir,max (km S I ) 


0.739±0.22 
[0.882±0.26] 


0.711±0.21 
[0.811±0.24] 


0.441±0.13 
(0.287±0.085) 


0.462±0.14 
(0.300±0.089) 


0.462±0.14 


0.352±0.10 




f7 N 2 D+ /^c^ir.max 


0.554±0.17 


0.573±0.17 


0.565±0.17 


0.902±0.28 


0.736±0.22 


0.853±0.26 


0.697±0.088 




[0.464±0.14] 


[0.502±0.15] 


(0.871±0.27) 


(1.39±0.44) 






[0.670±0.086] 
([0.803±0.11]) 


■Rc.vir.max (0.01 pc) 


7.74±2.6 
[8.16±2.8] 


8.44±2.9 
[8.81±3.0] 


4.33±1.5 
(1.82±0.62) 


3.18±1.1 
(1.34±0.46) 


4.81±1.6 


3.08±1.0 




fic/Kc,vir,max 


1.06±0.42 
[1.00±0.40] 


1.04±0.41 
[0.994±0.39] 


1.03±0.41 


1.05±0.41 


1.07±0.42 


1.05±0.42 


1.05±0.17 
[1.03±0.17] 


£ c ,mm (gmr 2 ) (30%) 


0.029 
[0.095] 


0.037 
[0.082] 


0.031 


0.031 


0.011 


0.018 




M Cjmin (Mq) (50%) 


2.92 
[9.56] 


4.26 
[9.44] 


0.926 
(0.165) 


0.522 
(0.0927) 


0.44 


0.286 




n H ,c,min (10 s cm- J ) (36%) 


0.368 
[1.21] 


0.439 
[0.973] 


0.722 
(1.71) 


0.962 
(2.28) 


0.221 


0.575 




^c^vii^min (km S 1 ) 


0.391±0.12 
[0.567±0.17] 


0.413±0.12 
[0.532±0.16] 


0.262±0.078 
(0.170±0.050) 


0.251±0.074 
(0.163±0.048) 


0.217±0.064 


0.190±0.056 




(T N 2 D+ /°"c,vir,min 


1.05±0.32 


0.986±0.30 


0.951±0.29 


1.66±0.52 


1.57±0.47 


1.58±0.48 


1.30±0.17 




[0.721±0.22] 


[0.765±0.23] 










[1.21±0.16 






(1.46±0.45) 


(2.56±0.80) 






([1.44±0.20 


) 


-Rc,vir,min (0.01 pc) 


2.17±0.74 


2.84±0.97 


1.53±0.52 


0.940±0.32 


1.06±0.36 


0.899±0.31 




[3.38±1.1] 


[3.79±1.3] 


(0.645±0.22) 


(0.396±0.13) 








Rc/ -RcjVir.min 


3.77±1.5 


3.08±1.2 


2.91±1.1 


3.56±1.4 


4.86±1.9 


3.61±1.4 


3.63±0.59 



[2.42±0.95] [2.31±0.91] 



[3.28±0.54] 



Si. 34mm (mjy) 

Si.34mm/^ (MJy/sr) 

Sc.mm (g cm" 2 ) 



6.94±0.72 
8.25±0.86 

-,0.321 



26.7±0.77 
27.7±0.80 
0.542, 1 ' ,JS 



5.08±0.63 
11.1±1.4 



0.218» ;tg 



3.70±0.30 
I I. ; : 1.2 

(°- 842 nl S ) 
8.74#§ 
(20.7| 1 ;D 



3.08±0.48 
3.10±0.48 

n 40b. U3 
z '^ z 0.992 



1.05±0.31 

2.66±0.79 

TTTTTUB 



0.161 



m 



M 22 

20.8 



0.0521 



m 



ML. 



Mc, mm (Mq) 



16.2: 



6.83 



62.5. 



nH,c,mm (10 b cm 3 ) 

ff c,vir,mm (km S ^) 



4.12 
1.10 



6.51^ 



0.826 



0.296 



5 07 W 
2-69 

0.427«|| 

(a27 i:l ) 



1 67 ;i - 41 
i - D, 0.766 



2.05 



1 oo2. 4b 
1 •■"<). 635 



0-60« 
[0-648»-l] 



0.808A;ii 3 



mrrrr 



U.OOOg 231 



0.248^;^ 



0.436, 

(0 ' 283 eP ) 



T^IS — 

0.736 



T777 — 

0.865 



-7X826^1- 
[0.813°;fp] 



%D+/ ff - 



c,vir ,mm 



0.681 
[0. 



0.487 



0-5048.364 

[0.4778 : ||1] 



0.584^ 



0.956 



1.02, 



1.21 



[4.4 ;gi 

1 fin 2 * ' ;) 



0-804A:^ 8 
[0-897A-fej,] 



(0-899A|2i) 

A nfiS-42 



^S 1 



9 /iQ^.yb 



^c,vir,mm (0.01 pc) 



4.06S.55 

(1.71?;^) 
j ijjTTra — 



2.83f-l 
sttri 



1.53,1 



z - u 'l.20 



1.22 



1.48H» 



C,vir,mm 



1.18A 



2.12' 



a Properties derived from BTK MIR extinction map for Cl-N & S, including sample averages, are shown inside "[...]". 

b Properties derived assuming the 1.56 kpc distance estimate to IRDC F (Kurayama ct al. 2011), including sample averages, are shown inside 
"(•••)"• 
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the core). Thus in the case of Cl-S, where cr N2D +/ cr c,vir,mm = 0.477 (i.e. with <j>B = 2.8), virial equilibrium would require 
(</>B,vir/2.8)~ 3 / 8 = 0.477, i.e. (j>B,vir — 20.1, corresponding to mA, V ir = 0.282 and -B c ,vir = 975 fiG. Analogous results for all 
the cores are shown in Table [3] For the cores that appear to be sub- virial in the fiducial analysis, a moderately stronger 
background magnetic field (by factors of a few) is required for virial equilibrium. Note that the ratio of i? c /i? c ,vir is only 
weakly affected by a change in ips, but stronger large-scale magnetic field support should lead to increased flattening of 
the core along the direction parallel to the field lines. 

Could the masses of the observed cores be set by their magnetic field strengths? For an ellipsoidal core of length 2Z C 
along the symmetry axis and radius R c normal to the axis, the magnetic critical mass, i.e. the maximum mass that can 
be supported by a magnetic field, -B c ,crit, is (Bertoldi & McKee 1992) 




If we equate M c b = M c mm for our cores (assuming R c = Z c ), we can estimate the values of -B c ,crit that apply to their 
particular masses and densities. We find that these values are similar to £> c .vir for all six cores, even though there is a 
range of required Alfven Mach numbers, ?T£A,vir 

of a factor of seven. 

Unfortunately, it is difficult to observationally determine magnetic field strengths in IRDC cores and we do not have any 
direct constraints for these particular sources. Measurement of Zeeman splitting from molecules such as CN (Falgarone 
et al. 2008) require relatively strong lines, whereas the observed emission is typically quite weak. From 14 regions, with 
average density of rin = 9 x 10 5 cm -3 , Falgarone et al. derived a median value of the total field strength of 560 fiG. 
Such values are similar to those needed for our N2D + cores to be in virial equilibrium. Our most massive core requires a 
moderately higher value (and at a slightly lower density), but since massive starless cores (that will form massive stars) 
are rare objects, it is quite possible they require relatively unusual (stronger) magnetic field strengths. 

Magnetic field morphology can be measured via dust continuum polarization, but there are large uncertainties in 
estimating field strengths via the Chadrasehkhar-Fermi method. Relatively order magnetic field morphologies, and thus 
relatively strong fields strengths, have been claimed to be present around some massive protostars (e.g. Girart et al. 
2009), which is additional indirect evidence that dynamically important fields are also present at the earlier, starless core 
stage. 

3.5. Comparison to previous studies 

Csengeri et al. (2011) studied the dynamics of five massive protostellar cores (i.e. those that are already forming a 
protostar) in Cygnus-X, concluding that the velocity dispersions were smaller than the predicted MT03 values for a viri- 
alized core. They measured velocity dispersions from observations of H 13 CO + (1-0) (~ 4" angular resolution; 0.13 kms -1 
velocity resolution) and N2H + (l-0) (29" angular resolution; 0.13 kms" 1 velocity resolution). They assumed a single value 
of S G i = 1.7 gem -2 based on mm continuum emission. The main differences compared to our present study include: (1) 
they studied cores that are already forming stars; (2) they did not use deuterated species to trace the core gas; (3) they 
measured clump properties via mm continuum emission, rather than MIR extinction mapping. 

Pillai et al. (2011) studied the dynamics of cold cores in the vicinity of ultracompact HII regions, also concluding 
they were strongly sub-virial. They used NH2D emission (~ 4.5" angular resolution; 0.27 kms -1 velocity resolution) to 
measure velocity dispersion and 3.5 mm continuum emission to estimate mass. The main differences compared to our 



Table 3 

Dynamical Properties of N 2 D + Cores 3 



Core property (% error) 


Cl-N 


Cl-S 


Fl 


F2 


G2-N 


G2-S 



R c (0.01 pc) (20%) 

CT N,D+ (kms" 1 ) 
S cl (gem- 2 ) (30%) 

Afc.mm (Mq) 
"H,c,mm (10 5 cm" 3 ) 

tc.fi (10 5 yr) b 
B c (mG) (m A = 1) 
-R c , v i r , mm (0.01 pc) 
Rc/R c ,vir,mm 

°- c ,vir,mm (kms -1 ) 

cr N2D+ /^e,vix,mm 

</>B,vir 
mA,vir 
B C ,vir (mG) 
Bc,crit (pG) 



8.18 

0.409±0.031 
0.483 

16 - 2 fi 3 83 

2.05™ 
0.962iji 8 

15g223 
fi 4Q0'.§§3 

u.djz 452 

9.53" 
0.427; 

365' 



219 



84 2 5 61 



8.75 

0.407±0.019 
0.402 
62.5ifo 
6,43 g ^2 

97K390 
Z '°2Q2 

9.761 
0.897J 
0.854i 
0.477S 



20.1 
0.282! 



'0.345 

48.0 



975)i 67d 
736 



fiio 



4.46 
0.250±0.021 
0.217 
6-5li 3 7? 

5 - 07 2°§ 

O-eilglf 
0.427°9 e flg 



0.584! 

11.8; 
0.379; 

304 ; 

296' 



0.410 
20.0 



3.35 
0.417±0.044 
0.324 

4 74 9. 80 

8.74*™ 

0.465§ : 63| 

"HIP 74 
234 
9 qq.4.48 

Z.OO-r 70 

i - ±o 0,687 
n 4QK0.601 

A3 %w 



0.956, 
3.16 



0.075 
7.08 



1.899-1 

367™ 
009834 

ooz 234 



5.16 
0.340±0.018 
0.214 
2.42§-gi„ 
1.222 
1.25J 
94.41 
2.49? 
2.07 : 
0.333! 



3 8 



K7„! 



159 

L - uz ().73(i 
9 fi46.35 
z, "*0.992 

1-06^45 

89.3^1 

go 9I8O 



3.25 
0.300±0.020 
0.194 

0.826q;o9 6 

•3.41 



1.673 
1.07, 

92.0, 
1.53, 
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present study include: (1) Many of their cores are unresolved, with estimated deconvolved sizes that are much smaller 
than the angular resolution of their observations - yet all mm emission, i.e. mass, is assumed to originate from inside 
these effective radii. On the other hand, all of our cores are resolved. (2) No allowance was made for surface pressure 
terms in the consideration of virial equilibrium. 

4. CONCLUSIONS 

We used ALMA Cycle observations to search for N2D + (3-2) emission from four IRDCs, selected to be dark at 
wavelengths as long as 70 /um. Strong detections were made in all cases, leading to identification of 6 N 2 D+ cores. We 
assessed their properties, including their surrounding clump envelopes, via MIR extinction maps. We also measured core 
masses via the detected 1.34 mm dust continuum emission. We regard this mass estimate as the most accurate for the 
cores (even though it still shows factor of 2 uncertainties), because of the difficulty of separating out core from clump 
material in the MIREX-derivcd mass surface density maps. The MIREX maps do however allow us to estimate the clump 
envelope properties, which are needed for the surface pressure term in the virial equation. 

We assessed the dynamical state of the cores by comparing to the Turbulent Core Model of MT03. For the sample 
of 6 cores the ratio of the observed to the predicted (virial equilibrium) velocity dispersion is 0.81o;fi, while the ratio 
of the observed to predicted size is 1.54}-2g. Thus the cores appear to be close to the predictions of the model, which 
assumes virial and pressure equilibrium and invokes magnetic fields such that the Alfven Mach number is uia = 1. There 
is tentative evidence that the cores are slightly sub-virial compared to the fiducial model, especially in the case of the 
more massive cores Cl-S and Cl-N, with <~ 63 M Q and <~ 16 M , respectively. However, these cores could be close to 
virial equilibrium if they are threaded by moderately stronger large-scale background magnetic fields with strengths up 
to ~ 1 mG, implying that rriA — 0.3 — 0.4. 

To prevent fragmentation of the cores requires similar field strengths, which we have evaluated by equating observed 
core mass with the magnetic critical mass given its observed density. This may indicate that magnetic fields are dynam- 
ically important in setting the mass function of cores, as has been suggested by Kunz & Mouschovias (2009), especially 
preventing the fragmentation of the massive cores that eventually form massive stars. Magnetic suppression of fragmen- 
tation appears to be the most likely mechanism operating in IRDCs, since the cold temperatures (~ 10 K) indicate that 
strong radiative heating from surrounding protostars to raise the thermal Jeans mass (Krumholz & McKee 2008) is not 
occurring. Magnetically mediated massive star formation would not require a minimum mass surface density threshold, 
such as the £ c i ~ 1 g cm~ 2 limit proposed by Krumholz & McKee. Indeed the clump medium immediately surrounding 
Cl-N and Cl-S only has S c i ~ 0.4 g cm -2 and the clumps in the IRDC sample of BT12 tend to have even lower values. 

This pilot study with ALMA revealed two relatively massive starless cores out of a total sample of six detected objects. 
Further observations of other regions are needed to build a larger statistical sample. Observations to constrain the 
magnetic field strengths in these objects are also needed. 
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acknowledges support from NSF CAREER grant AST-0645412; NASA Astrophysics Theory and Fundamental Physics 
grant ATP09-0094; NASA Astrophysics Data Analysis Program ADAP10-0110. SK acknowledges support from an NRAO- 
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